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We study the properties of spectra and eigenfunctions for a chain of 1/2— spins (qubits) in an 
external time-dependent magnetic field, and under the conditions of non-selective excitation (when 
the amplitude of the magnetic field is large). This model is known as a possible candidate for 
experimental realization of quantum computation. We present the theory for finding derealization 
transition and show that for the interaction between nearest qubits, the transition is very different 
from that to quantum chaos. We explain this phenomena by showing that in the considered region 
of parameters our model is close to an integrable one. According to a general opinion, the threshold 
for the onset of quantum chaos due to the interqubit interaction decreases with an increase of 
the number of qubits. Contrary to this expectation, for a magnetic field with constant gradient 
we have found that chaos border does not depend on the number of qubits. We give analytical 
estimates which explain this effect, together with numerical data supporting our analysis. Random 
models with long-range interactions are studied as well. In particular, we show that in this case the 
derealization and quantum chaos borders coincide. 



PACS numbers: 05.45Pq, 05.45Mt, 03.67Lx 



I. INTRODUCTION 

In recent years much attention has been paid to the 
idea of quantum computation [jjj. The burst of inter- 
est to this subject (see, for example, [^]-|| and references 
therein) is caused by the discovery of fast quantum al- 
gorithms for the factorization of integers pi and for the 
effective searching of items in a database |3|,|7| . These al- 
gorithms demonstrate the effectiveness of quantum com- 
puters in comparison with the classical ones. Nowadays, 
there are different projects for the experimental real- 
ization of quantum computers, as well as experimental 
results with few-qubit systems (see j8|) and references 
therein) . 

Main theoretical suggestions for the experimental im- 
plementation of the quantum computation are based on 
interacting two-level systems (qubits ). It is clear that 
one of the most important problems from the viewpoint 
of the stability of quantum operations, is a destructive 
role of different kinds of errors. In the first line, one 
should refer to finite temperature effects and the interac- 
tion with an environment |J. However, even in the case 
when these effects can be neglected, there are dynamical 
effects of the interqubit interaction, which may influence 
a quantum computation. On one hand, the interaction 
between qubits is necessary for the realization of quan- 
tum computation, on the other hand, it may result in a 
kind of destruction of the coherence in the evolution of a 
system. 

The latter subject of the dynamical decoherence is di- 



rectly related to the so-called quantum chaos which is 
nowadays widely discussed in application to atoms, nu- 
clei, quantum dots and other physical systems (see, for 
example, |I(J and references therein) . One of the latest 
developments in the theory of quantum chaos refers to 
the interaction between Fermi-particles in isolated sys- 
tems. The core of this approach is the perturbation the- 
ory for many-body states, which takes into account a two- 
body nature of the interaction. Specifically, it was shown 
|]llf that if the two-body random interaction between par- 
ticles exceeds some critical value, fast transition to chaos 
occurs in the Hilbcrt space of many-particle states (see 
also and reviews ||@). 

In dynamical systems such as complex atoms |18[] , mul- 
ticharged ions JlSfl , nuclei |2(J and spin systems [^1 22 
quantum chaos gives rise to a very complicated structure 
of highly excited states, and to specific correlations in 
the energy spectra, described by Random Matrix The- 
ory (RMT) (see, for example, Jl6|]). As a result, closed 
dynamical systems with relatively small number of in- 
teracting particles can be well described by a statistical 
approach, see discussion and references in p3[ . 

Recently, quantum chaos theory has been applied to 
a simple model of quantum computer p4| ] chosen in the 
form of L interacting qubits. Numerical data have shown 
that for a strong enough interaction between qubits the 
onset of quantum chaos is unavoidable. Although for 
L = 14 — 16 the critical value J cr for the quantum chaos 
border was found to be quite large, with an increase of 
L the border decreases as J cr ~ 1/L p5|,p4j. From the 
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viewpoint of the standard approach for closed systems of 
interacting particles, the decrease of chaos border with an 
increase of qubits looks generic. This poses the question 
of the relevance of quantum chaos to quantum computa- 
tion 

In our recent paper [ p9[ we have studied the errors 
which appear in the evolution of ID Ising nuclear spins in 
rotating magnetic field. This model was suggested for an 
experimental realization of a quantum computer p(i| , |3l| . 
The main attention in |29[ has been paid to the region 
of parameters, most suitable for the preparation of an 
initial many-body state needed for further application of 
quantum protocol (sequence of time-dependent magnetic 
pulses in a prescribed algorithm of quantum computa- 
tion). It was shown that even for a very large interqubit 
interaction, the errors turn out to be very small, thus 
demonstrating that the influence of quantum chaos can 
be neglected. 

An analysis of the stationary Hamiltonian describing 
the system during a single magnetic pulse has been per- 
formed in Q . Specifically, the general approach of quan- 
tum chaos theory has been applied, in order to under- 
stand the conditions for the onset of quantum chaos. The 
model we considered assumed that qubits (nuclear spins) 
are placed in a strong magnetic field with constant gra- 
dient along the direction of the spin chain. The gradi- 
ent of the magnetic field provides a "labeling" of qubits. 
Namely, each spins has different Larmor frequency, ujk- 
This allows one to provide a selective addressing to each 
qubit by applying resonant rf pulses. The main interest 
was in the influence of the magnetic field on the proper- 
ties of eigenstates and energy spectra. It was unexpect- 
edly found that the constant gradient magnetic field gives 
rise to the independence of the critical value J cr on the 



number L of interacting c 
ena has been explained in 



:ubits. This striking phenom- 
[ 2j| analytically and confirmed 
numerically, thus giving a new insight to the problem of 
quantum chaos in the models of quantum computers. 

In this paper we present the full theory which explains 
the properties of energy spectra and many-body states 
of the model of Ref . |2J|] , together with numerical data 
obtained in a broad region of the model parameters. The 
structure of the paper is as follows. In the next section 
we describe the model, discuss the region of parameters 
of our interest, and briefly analyze the structure of the 
Hamiltonian matrix in the z-representation. In Sect. Ill 
we study global properties of the energy spectrum, pay- 
ing main attention to the band structure of the spectrum 
and to the level spacing distribution P(s) for the central 
energy band. 

Section IV is the core of the paper, here both the delo- 
calization border and the condition for the onset of quan- 
tum chaos are studied. The consideration has been made 
by making use of the mean-field representation which is 
very convenient from the theoretical viewpoint. One of 
two main goals of this section is that these two borders 



are very different in the model with nearest interaction 
between qubits. Another important result is that the de- 
localization border turns out to be independent of the 
number of qubits for a gradient magnetic field. Theoret- 
ical estimates obtained in this Section serve as a guiding 
line to treat all numerical data. 

In Section V we investigate numerically the structure 
of eigenstates in the ^-representation, by relating the 
data with the theoretical predictions. Section VI is de- 
voted to some modifications of the model, namely, we 
analyze the influence of randomness in the interqubit in- 
teraction. Our main question is how statistical properties 
of the system depend on the range of the interaction be- 
tween qubits. Specifically, we study random interaction 
between all qubits (A -interaction), as well as between 
four nearest qubits, by comparing the results with those 
obtained for the model with the interaction between two 
nearest qubits ( A-interaction) . 

General discussion is presented in the last Section VII. 
One of the problems we discuss here, is the concept of 
the quasi-integrability of our model for the A-interaction. 
We show that for the region of parameters of our interest, 
the model is close to the integrable one. This explains 
why the delocalization and chaos borders do not coincide 
for the A-interaction. We also analyze the role of the 
magnetic field. In particular, we give analytical estimates 
which show that for the homogeneous magnetic field the 
delocalization border has generic L-dependence discussed 
in |24|]. On the other hand, for the magnetic field with 
an increasing gradient, analytical estimates predict that 
the delocalization border increases with an increase of 
the number of qubits. 



II. THE MODEL 

The model describes a 1-dimensional chain of L inter- 
acting distinguishable 1/2-spins in an external magnetic 
field. Schematically, these spins (qubits) can be repre- 
sented as follows, 
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Here B z stands for a constant part of magnetic field ori- 
ented in the positive z-direction, and each qubit occupies 
one of two single-particle states with the energy 1/2 (po- 
sition "up") or —1/2 (position "down"). One can see 
that the total number A of many-body states which are 
generated by this chain (quantum register), is A = 2 L . 

The dynamics of this model (quantum computer pro- 
tocol) is due to a sum of p — 1,...,P time-dependent 
rectangular pulses of a circular polarized magnetic field 
rotating in the x, y-plane. Each of the pulses has own 
amplitude 6^, frequency v p , phase tp p , and lasts during 
the period T p — t p+ i — t p . Therefore, the total magnetic 
field during one pulse can be written as follows, p9||, 
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B(t) = [b\cos{v p t + <p p ),-b\ sm{v p t + <p p ),B z ], (1) 
The Hamiltonian of this system has the form, 

n = - eW| + 2 E Animh 

k=0 n>k 

\ E e P (t)n p E ( e-^-^ei- + e «v+^ P/ + ] , 



p=l 



(2) 



where the "pulse function" Q P (t) equals 1 only during 
the p-th pulse of the length T p . The quantities Jk, n 
stand for the Ising interaction between two qubits , uj k are 
the frequencies of spin's precession in the B z — magnetic 
field, fl p is the Rabi frequency corresponding to the p- 
th pulse. The operators Ij~ arc defined by the relations 
if = I x k ± il v k , and K' V ' Z = (l/2)<' a,z , the latter being 
the Pauli matrices. 

Below we consider the properties of the system during 
a single p-th pulse. The corresponding Hamiltonian can 
be written in the coordinate system which rotates around 
z-axes with the frequency v p . Thus, for the p-th pulse, 
our model can be reduced to the stationary Hamiltonian, 



= - L EIK 

k=0 

2 E Jk^I^n 
n>k 



v p )ll + tt p (cos(fi p I£ - sin Lp p I\ 



I z ] 



(3) 



which describes the evolution of the model for t p < t < 
tp+i- 

The regime of quantum computation corresponds to 
the following range of parameters: fl p <C J kl n *C Su k 
u>k, where 5uj k — — ui k \ p9| (the so-called selective 

excitation) . In this regime, each pulse acts selectively on 
a chosen qubit exciting a resonant transition. The in- 
equality, Q p <C Jk,n, provides a separation between reso- 
nant and non-resonant transitions for the same selected 
qubit. The inequality, J k , n ^ SiOk, provides a separation 
of transitions for a given qubit from the transitions for 
neighboring qubits. In this paper we consider another 
important regime of non-selective excitation which is de- 
fined by the conditions, tt p 3> 5uj k 3> J, see details in 
p9[ . This inequality provides the simplest way to pre- 
pare a homogeneous superposition of 2 L states needed 
for implementation of both Shor and Grover algorithms. 

In what follows we assume, for simplicity, tp p = 7r/2, 



and put Op 



fl and Vr, — 



Our main interest is in 



the nearest neighbor interaction (N- interaction) between 
qubits for two different cases, the dynamical one when all 
coupling elements are the same, J k , n = J <5n,fc+i, and the 
case when all values Jk,k+i are random (random model). 
However, we will also analyze other cases with differ- 
ent kinds of interaction and compare results with those 



for the iV-interaction. In contrast to the previously dis- 
cussed model |Q with homogeneous magnetic field, be- 
low we consider the magnetic field which depends on the 
position of the fc-th qubit. Therefore, we assume that 
the spin frequencies uj k are slightly dependent on k (with 
8w k < Wfc). 

For the dynamical N— interaction, the Hamiltonian (|^) 
takes the form, 



L-l 



L-2 



H = 



(4) 



fe=0 



fc=0 



where S k — uj k — v. In the ^-representation the Hamil- 
tonian matrix of size N = 2 L is diagonal for O = 0. For 
O ^ the off-diagonal matrix elements are H k ^ n = iO/2 
for n > k, and H nk = H^ n . When calculating the matrix 
elements of the Hamiltonian (|j) we have used the stan- 
dard rules in order to find the action of the operators I| 
and I\ on the states |fc) and 1 77,) , 

i?\...o k ...) = h...o k ..,), ii\...i k ...) = -h..i k ...), 



- l -\...i k ...), q\...x k ...) = - l -\ 

The matrix turns out to be very sparse, and it has 
quite specific structure (see Fig.l) in the basis which 
is reordered according to an increase of the number s 
written in the binary representation, s — Il—i, H) *o 
(with i s = or 1, depending on whether a single-particle 
state of i— th qubits is the ground state or the excited 
state). In what follows, we call this representation as the 
z — representation. 



In> 
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FIG. 1. Structure of the Hamiltonian (^) for L = 8 in 
the z -representation. Dots stand for the matrix elements 
different from zero. 



III. GLOBAL PROPERTIES OF THE ENERGY 
SPECTRUM 

For the further analysis, it is important to understand 
the global structure of the energy spectrum. In what fol- 
lows, we concentrate our attention to the case when the 
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magnetic field has a constant gradient along the chain of 
qubits, Wk = Wq + ak with a > 0. Other cases will be 
briefly discussed in Sect. VII. 

A. Band structure 

Without the interaction between qubits, J = 0, the 
energy spectrum of the model (Q) consists of L + 1 bands 
of finite width for a ^ 0, separated by big gaps of size 
f2 ^ uik ■ In Ref . |^j],^8| it was numerically found that the 
width AE(il, J = 0) of the central band decreases with 
an increase of £1 as A^/fl. Our analytical estimates show 
that for L even, the bandwidth is given by the relation 
(A.E)i = L 2 a 2 (L - l)/80, (see details in Sect. IV). This 
dependence also occurs for a relatively weak interaction 
J 7^ 0. However, when the interaction exceeds some crit- 
ical value J s , the band widths turn out to be practically 
independent on 57, see the data for the central band in 
Fig. 1 




FIG. 2. Dependence of the central bandwidth AE on f2 
for different values of J. The data are shown for L — 10, 
ujk = u>o + k, v = u)q and a = 1. The full straight line is 
the theoretical expression for J = 0. The horizontal lines on 
the right hand side of the figure correspond to the analytical 
expression for (AE)2 for the case of J > J s , (see in the text). 

The bandwidth (AE)2 for the interaction strength J 
larger than the critical value J s can be also estimated an- 
alytically as {AE) 2 = (L- 2)aJ (see Sect. IV). The cor- 
respondence between the analytical estimate and numer- 
ical data was found to be quite good. If the band-width 
AE is larger than (AE)2, the first expression (A_E)i for 
AE(fl) dominates. Contrary, if the band-width (AE)2 
defined by the interaction J is larger, it determines the 
actual band-width AE which is independent of Q. 

One should stress that the above consideration is valid 
for the case when the bands are not overlapped. One can 
expect that for sufficiently strong interaction between the 
qubits, the band structure of the energy spectrum dis- 
appears. The overlapping of the central band with two 
other bands is shown in Fig.0 where the edges of the cen- 



tral and the nearest bands are plotted in dependence on 
the interaction J for the fixed value fl — 100. One can 
see that for J > Jb ~ 15 the bands are overlapped, there- 
fore, a change in the properties of the system is naturally 
expected. The critical value J& for the overlapping of the 
bands is estimated in Sect. IV as well. 

100 ^ 



E 




-100 r 

10 20 30 

J 

FIG. 3. Energy bands as a function of J. Only the 
central band (full line) and its neighbors (dashed lines) are 
shown, thus demonstrating the band overlapping for a rela- 
tively large interaction. The parameters are L = 10, fi = 100, 
ujk — + k, v — luq and a = 1. 



B. Level spacing distribution 

Let us now analyze the distribution P(s) of spacings 
s between nearest-neighbor energy levels inside the cen- 
tral energy band (note that s should be normalized to the 
mean spacing between levels). This quantity is often used 
in the theory of Quantum Chaos as a detector of chaotic 
properties of a system. Specifically, for systems with reg- 
ular motion in the classical limit, the distribution P(s) is 
generically close to the Poisson (apart from 1-dimensional 
systems where P(s) is highly non-generic and can be of 
any form) . In the other limit case of a completely chaotic 
motion, in the corresponding quantum systems the distri- 
bution P(s) has the so-called Wigner-Dyson (WD) form 
which is characterized by the level repulsion for small 
spacings, s <§; 1 (P(s) ~ s, s 2 , s 4 , depending on the sym- 
metry of a system, see, e.g. (lo|). 

Numerical data for P(s) for different values of the in- 
teraction strength J, summarized in Fig^, manifest the 
transition to the WD-distribution. Note that for small 
values of J (see Fig.[|a-c) the distribution P(s) reveals a 
strong deviation from the Poisson. Specifically, one can 
detect a clustering of energy levels for very small s, that 
results in a huge peak in the distribution at the origin 
s = 0. The presence of this peak indicates that for weak 
interaction our model is highly non-generic and may be 
compared to integrable ID-models. With an increase of 
the interaction, data for P(s) reveal, first, a transition 
to the Poisson, and after, to the WD-distribution. The 
more detailed analysis of the data shows that the tran- 
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sition from Poisson to WD-distribution occurs when the 
central energy band starts to overlap with the nearest 
bands. 



P(s) 



P(s) 




this way the "mean field " Hq describes a regular part of 
the total Hamiltonian, and the term V, describing the in- 
teraction between the particles, is responsible for chaotic 
properties (if any) of the system. Such a mean-field ap- 
proach is typical in the study of chaotic properties of 
complex atoms and heavy nuclei [|18|]2C 



A. Derealization border 

Now, one needs to represent the Hamiltonian (|^) in 
the basis in which it is diagonal in the absence of the in- 
teraction (J = 0). In this representation (corresponding 
to the rotating basis) the Hamiltonian Hq can be writ- 
ten as a sum of L individual Hamiltonians H k describing 
non-interacting quasi-particles |28|, 



FIG. 4. Level spacing distribution P(s) in depen- 
dence on the interaction J = 0,0.0002,0.1,1,10,100 de- 
noted by (a, b, c, d, e, f) respectively. Other parameters are 
L = 12, n = 100 , uj k = + k, k = 1,...,L, uj = 100, 
v = uo and a = 1. Full curve is the Poisson, dashed curve 
corresponds to the Wigner-Dyson distribution. Eigenvalues 
are taken from the central energy band only. 



IV. THEORY 

Let us now discuss our model ^ from the viewpoint of 
the standard approach to interacting particles in isolated 
systems (see, for example, |lj]]2^] and references therein). 
In this approach the Hamiltonian is written in the form 
H = Hq + Vq where Hq stands for non-interacting par- 
ticles, and Vq describes a two-body interaction between 
particles. The onset of chaos is usually meant to occur 
when the strength of the interaction Vq exceeds the mean 
energy spacing df between those many-particle states 
which are directly coupled by the interaction. It is im- 
portant to note that this spacing is much larger than the 
mean level spacing D between many-particle states. In- 
deed, while the total density p = D^ 1 of states increases 
exponentially with the total energy, the density pf = dj 1 
increases only algebraically (for details see, e.g., [19]). 

In order to apply this approach to our model (4|), one 
needs, first, to present the Hamiltonian as a sum of the 
"unperturbed" Hamiltonian Hq, and the part V = JVq 
that stands for the interaction between particles. In our 
case the Hamiltonian (0) can be rewritten in the form, 



H = H + JVq. 



(5) 



where 



L-l 



L-2 



Ski 



JMfc 



k=0 



Vq 



-2£4*4V (6) 



fc=0 



As one can see, the Hamiltonian Hq stands for a kind 
of mean field which absorbs the Sl-dependent term. In 



L-l 

Ho = £ H k = 

k=0 



L-l 

E 

k=0 



S 2 k + Wl z k . (7) 



Correspondingly, in the basis of Hq the interaction Vq be- 
tween quasi-particles has the form, Vq — Vdiag + Vband + 
V ff , where 

Vdiag = -2J2hh+iIk^k+i > 



V, 



band 



- -^J2 a kak+iitik+i ; 



Voff = 2 E {a k bk + il y k l z k+1 + a k+1 b k r k I v k+1 ) 

k 



where 



n 



«k 



(8) 



(9) 



From Eq.(@) one can see that the energies e k of quasi- 
particles (or, the same, energies of single-particle states 
determined by the Hamiltonian Hq) are given by the ex- 
pression, 



± 



(10) 



Note that this relation is valid for any kind of magnetic 
field B z (any dependence 8k), including the homogeneous 
magnetic field (S k = constant). 

Let us now consider the constant gradient magnetic 
field (6k = ak) for large values of f2 ^> S k . In this case, 
one can write an approximate relation for e k , 



a- 



±- 



1 



-'A: 2 



2il 



(11) 



This expression allows one to find global properties of 
the unperturbed (J = 0) energy spectrum, briefly dis- 
cussed in the previous section. Indeed, for large values 
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of £1 (more correctly, for fl ^> ak) one can see that the 
spectrum has a band structure, with the bands centered 
at 0, ±ft, ±2fi,...,±iO. 

The central band is defined by such locations of quasi- 
particles in the single-particle spectra defined by e^, for 
which an equal number L/2 of quasi-particles have pos- 
itive and negative values of (for an even number L of 
qubits). Therefore, the total number N c b of many-body 
states in the central band is given by the total number of 
combinations of N objects having half positive and half 
negative values, 



(L/2)!(L/2)! 



(12) 



One can also see that for 8k = ak and J = 0, the size of 
the central energy band is given by twice the maximum 
energy inside the band, 



{AE) cb = 2E { ™ ax ^ = 2 



' L-l L/2-1 

E k2 - E k2 

k=L/2 k=0 



(13) 



_ L 2 (L - l)a 2 

Now, let us discuss the structure of the Hamiltonian 
matrix determined by the off-diagonal terms (|J). One 
can see that in the unperturbed basis the term Vdiag 
is clearly diagonal. The action of Vband, is much more 
complicated. Let us consider, for simplicity, the central 
band. Each operator It flips the k-th spin. Since the 
interaction is two-body, we should consider the action of 
1 u P on states as |...,0 fe+ i, l fc , ...) , |..., l fc+ i, fc , ...) , 
|...,0 fe+ i,0fc, ...) ,|..., lfe+i, lfc, ...). First two kinds of the 
states upon the action of Vband still remains in the same 
central band since the number of 0's and l's is conserved. 
Second pair of states increases (or decreases respectively) 
the number of l's of two units: that means that such a 
coupling refers to a next to nearest energy band (near- 
est bands differ by plus/minus one l's). As a result, 
one can conclude that the term Vband stands for the in- 
teraction both inside the central band, and between the 
next-neighbor energy bands. 

In the same way it is easy to understand that the term 
V ff give rise only to the off-band interaction, to be more 
precise, to a coupling between nearest bands. The struc- 
ture of the Hamiltonian in the mean-field basis is shown 
in Fig.| 

For a relatively weak interaction, the eigenstates in 
the mean-field basis defined by the unperturbed Hamil- 
tonian Hq, are delta-like functions with an admixture of 
other components with small amplitudes. In this case 
one can speak about the localization of eigenstates in the 
unperturbed basis. With an increase of the interaction 
strength, the number N pc of basis components with large 



amplitudes (number of principal components) increases. 
According to the theory of interacting particles, the tran- 
sition from strongly localized (N pc « 1) to delocalized (or, 
extended) states (with N pc 3> 1) occurs very fast with an 
increase of the interparticle interaction. For this rea- 
son, one speaks about the delocalization transition (in 
the finite-size basis), see, e.g. pTl and references therein. 
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FIG. 5. Structure of the Hamiltonian matrix in the mean 
field basis for the iV-interaction; black points stand for matrix 
elements whose modulus is larger that 10 -6 . Here is L — 8, 
Q = 100, J = l,u k =u + k. 

Generically, in the models with two-body random in- 
teraction V between particles pq |, extended eigenstates 
with large N pc turn out to be chaotic. By this term we 
mean the situation when the components of the extended 
states can be treated as random and independent quan- 
tities. Similar situation (the onset of quantum chaos) oc- 
curs in many dynamical systems with complex enough in- 
teraction, such as many-electron atoms and heavy nuclei 
Pq , P0| . In these systems, the delocalization transition 
coincides with the transition to chaos, and is determined 
by the condition V > df (with V as a typical interaction 
strength, and df as the mean energy distance between 
directly coupled many-body states). 

Let us now discuss the delocalization transition in our 
model, keeping in mind that it can be different from the 
transition to chaos. As it will be shown, our model with 
the iV-interaction manifests a quite unexpected phenom- 
ena, namely, the above two transitions turn out to be 
very different. 

We start with the estimate of the mean level spacing df 
in the central energy band of our model (|^) , between the 
many-body states coupled by the interaction (||). The 
energy spacing df can be estimated as the ratio 



(14) 



where Mf is the number of many-body states coupled by 
Vband inside the energy interval (AE)f. In fact, Mf is the 
mean number per line of non-zero off-diagonal elements 
in the total Hamiltonian (gj). 



G 



In order to estimate Mf, we note that the interaction 
Vband in the central band can only couple those many- 
body states having an equal number (L/2) of spins "up" 
and "down" (for an even number L of qubits) . The min- 
imal value of Mf = 1 corresponds to the state 

|0 L _x,..,0,l,...,l ), 



and the maximal one, Mf 
state 



L — 1, corresponds to the 



>L-l, J-t-2 



) 0l_3, 1l_4, Oi, 1 ). 



AE f 




L 

FIG. 6. Numerical calculation of (AE) / for both N- (open 
circles) and A-interaction (full circles), see Sect. VI) showing 
the //-dependence as found analytically. Here is a — 1. 



Indeed, in the first case there is only one possibility to 
change to 1 , and 1 to for the nearest qubits. And 
in the second case, there are L—\ such changes, each of 
them corresponding to the nearest neighbor interaction 
with no change in the total number of spins "up" and 
"down" . Therefore, one can estimate the average value 
Mf as Mf « L/2 which is in a very good agreement with 
direct numerical check. 

One should stress that the energy range (AE) t within 
which the many-body states are coupled, is much less 
than the total energy width (AE) c b of the central band 
determined by Eq.(§J|). The value of (AE) f can be 
estimated as the maximal difference between energies 



E {2) = J2[ 2) e k and E, 
states \ipi) and |?/>2) of Hq, having the matrix clement 
(ipi\Vband\ip2) different from zero. If we consider only the 
coupling inside the central band we can find these values 
E^ and E^ by observing that the maximal energy is 
obtained by flipping the outermost spins. Application of 
to tne state 



(i) _ 



£fc of two many- body 



l^i 



gives rise to the state 



|U-i,0 L _ 2 ,. 



As a result, for L > 1 we have, 



d f = 



(A£)j 
M f 



2a^_ 



(15) 



The mean spacing df should be now compared with 
the typical value of the perturbation, V = JVq. The lat- 
ter can be found from Vband as V ~ J/2 (other terms are 
negligible for 3> 6k). Therefore, we finally obtain, 



(16) 



Surprisingly, the derealization border does not depend 
on the number of qubits, in contrast to the result of ]24| ] 
where J cr decreases as 1/L. The reason is the specific 
influence of a constant gradient of the magnetic field 
that results in the quadratic dependence e k oc k 2 for the 
single-particle levels of quasi-particles of the mean-field 
Hamiltonian H (see Eq.(fj])). 



Let us now compare the analytical estimate (16) with 
numerical data. The commonly used quantity to measure 
the number N pc of principal components in eigenstates, 
is the so-called inverse participation ratio, 



N pc (E) 



-i -i 



(17) 



IV2) = |0l_i,1l-2,...>; 

(L — 1 and L — 2 correspond to the states with the high- 
est values of single-particle energies £&). Thus, the energy 
difference E\^\ — E\^, 2 \ is given by 



Numerical results confirm this prediction very well, see 
Fig§ 



Here ip n (E) = (n\ip(E)) is the n — th component of a 
particular eigenfunction corresponding to the eigenvalue 
E. 

From Eq. ( p7| ) one can see that for equal values of the 
components of an eigenstate, ip n = 1/yN, the number 
of principal components is equal to the size of the basis, 
N pc = N. In another extreme limit of completely ex- 
tended and chaotic eigenstates, the value of N pc is equal 
to N/3. The factor 3 arises due to the Gaussian fluc- 
tuations of ipn which are generic in the case of strong 
quantum chaos (see, e.g. |l^] ). For localized states the 
value of N pc approximately gives the number of basis 
states effectively occupied by this eigenstate. 
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FIG. 7. The average number of principal components in 
the rotated basis for the eigenstates from the central energy 
band, as a function of J/ J cr for L — 8, 10, 12 and different 
values of f2. 

Numerical data for N pc computed in the mean field 
basis (where Hq is diagonal for J = 0, sec Eq.(H)) for 
the eigenstates taken from the central energy band, are 
given in Fig.J^ in dependence on J / J cr . It is clearly seen 
that below the delocalization border, J < J cr , there is 
a scaling dependence of N pc on L and 17 in accordance 



with the estimate (16). On the other side, for J ^s> J c , 



the value of N pc saturates to its maximal value iV c h/3 
in correspondence with random matrix predictions (here 
N c b is the total number of states inside the central en- 
ergy band, see (|l~2|)). The latter correspondence of the 
maximal value of N pc to N c b/3 is a strong evidence of the and 
quantum chaos in the model for a very large interaction. 



In the case of aL 3> J we are interested in, the unper- 
turbed (J = 0) many-body energy spectrum of Jis| ) is 
given by a sequence of degenerate levels separated by the 
spacing a. Due to a weak interaction J , each set of these 
degenerate levels spreads and creates the energy bands. 
In order to find the central energy bandwidth, we should 
consider the action of the interaction operator 

L-2 
k=0 

upon the states belonging to the central band The lat- 
ter assumption is an approximation: in z-representation, 
the central energy band can also contain few states with 
slightly different numbers of O's and l's. 

The action of each term in the interaction operator 
leaves the state as it is, multiplying it by a factor ± J/2 
depending on the presence of two close 11 and 00 or differ- 
ent 01 and 10. This results in the shift of the energy from 
its zero-value in the central band. Two configurations, 
I-) = |0,0,...,0,1,1,...,1) and |+) = |0,1,0,1, ...,0,1), 
should be considered which correspond to the maximal 
shift in the "negative" and "positive" directions. In this 
way we can safely say that such a bandwidth is given by 
the energy difference Ei + \ — Ei_\. It is easy to see that 
one has 



E\ 



\+) 



(L - 1) Jo 



(L - 3)Ja 



B. Chaos border 

In this section we study the transition to global chaos 
which is due to the overlapping of the energy bands. In 
order to obtain the condition for the band overlapping, 
one needs to find the band width and to compare it with 
the distance f2 between the bands. 

We have shown that in the absence of the interaction, 
the energy width of the central band can be estimated an- 
alytically, see Eq.(|l3|). Numerical data reported in FigJ|, 
show that with an increase of the interaction, the band 
width AE saturates to some value (AE) S which is in- 
dependent of fi. Therefore, we can estimate (AE) S by 
coming back to the z— representation of Hamiltonian (^ 
) where the Jl-depending term enters in the off-diagonal 
matrix elements only. By omitting this term, we can 
write the diagonal part, 



L-l 

£ 

fe=0 



L-2 

ofe/|+2J^J|Jf +1 

k=0 



(18) 



where the relations LUk 
taken into account. 



ujo + ak, and v = luq are directly 



thus giving 



AE = (L — 2)Ja 



for 



n > j 



(19) 



By equating two expressions Eq. (|l3|) and (|T^) we finds 
the transition point 



Jo 



L 2 a 

~87T 



(20) 



between the two dependencies for the band width AE. 

One can see that for J > Jo bands are overlapped if 
(L — 2) Jo > f2. That gives the critical value Jt, for the 
overlapping, 



(21) 



J b «— . 
aL 



subject to the condition J > J . By comparing Eqs.(|2£ 
and (|2l]), one gets the lower bound J > compati- 
ble with the above two constraints. 

On the other side, one can also have the band overlap- 
ping whenever J < J , if a 2 L 2 (L— 1) /80 > fl. Therefore, 
in this case the overlapping of the bands occurs for any J, 
if the number of qubits is large enough, L > 2(fi/a) 2 / 3 . 
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One should stress that overlapping of bands is not a 
sufficient condition in order to have the derealization 
of eigenstates. Indeed, the estimate Eq.([l6|) for the de- 
localization border J cr is derived for the central energy 
band only, therefore, it is not valid when bands are over- 
lapped. Therefore, one needs to start with the expression 
( |l4| ) and estimate (AE) f for the case when the energy 
spectrum is not a band-like. 

In order to do this, it is convenient to switch to the 
mean-field representation with the unperturbed Hamil- 
tonian Hq given by Eq.(^| ). The total size of the unper- 
turbed energy spectrum is now defined by the difference 
between the energies corresponding to the following lim- 
iting configurations, 

| T) = |1l-i,U-2,...,1i,1o) 

and 

I I) = lOi-i.O^.-.^Oo). 

But this is not what we need. Indeed, these two many- 
body states are not coupled by the two-body interaction 
(|J) . What we need to find, is the maximal energy change 
due to the action of the J-interaction. To do that we 
have to consider two states corresponding to the flipping 
of both the two uppermost spins, namely : 

| T) = |1l-i,U-2,...) 

and 

U) = |0i-i,0i-2,...)- 

The energy difference between such states is given by : 

(AE) j « 4 v / a 2 L 2 + Q 2 . 

Since the number of coupled states remains the same, 
Mf ~ L/2, one gets that in order to have the transition 
to delocalized states in the case of the overlapped bands, 
the typical value of the interaction has to be larger than 
d f = (AE)f/Mf, 

' ~ L/2 

or 

i a 

J> J c ~ — Va 2 L 2 + n 2 (22) 
Li 

One should notice that the two criteria (band over- 
lapping and transition to derealization), if both satis- 
fied, are expected to result in the onset of chaos. This 
conclusion is confirmed numerically, and is supported by 
analytical arguments. 

Indeed using data from Fig{| e.g L = 12, Q = 100, 
a = 1, one gets a chaos border J c w 130 well confirmed 
by the Wigner-Dyson distribution in Fig.fff. 



On the other hand, we have already seen that the sim- 
ple requirement to be in the delocalized regime (J > J cr ) 
without the overlapping of bands, does not give rise to 
chaos in our dynamical model with the nearest interac- 
tion. 

It is also easy to check that the conditions of the band 
overlapping for the case J < Jq = Q/8ax and L > 8x 
(with x = fl 2 /L 2 a 2 ), are not compatible with the der- 
ealization border J > J c — 16ay/l + x in the region of 
non-selective excitation, i> 1. This means that a rela- 
tively weak interaction does not lead to the derealization 
(and, therefore, to the chaos), in spite of the overlapping 
of the energy bands. 

V. STRUCTURE OF EIGENSTATES IN THE 
^-REPRESENTATION 

The analytical treatment we have performed in the pre- 
vious section, is based on the mean-field representation 
of our model, namely, when the Hamiltonian matrix is 
written in the basis of the "unperturbed" part Ho, see 
Eq.(§. This approach is natural for the theoretical study 
since the interaction is much less that the ^-dependent 
term (J <C f2), therefore, the interaction between qubits 
can be considered as a weak perturbation. 

However, the dynamical properties of the model are 
related to the ^-representation which is adequate to the 
experimental setup. For this reason we discuss below 
the structure of eigenstates of Hamiltonian (|4| ) in z- 
representation, in relation with the above analytical es- 
timates obtained in the mean-field approach. 

Since the most important question is about the role of 
the interqubit interaction, main attention is paid to the 
dependence of global properties of eigenstates on the in- 
teraction strength J. Typical structure of the eigenstates 
in the z-representation is shown in Fig.^ for different val- 
ues of J. First, one should note that in this basis all 
components of eigenstates in the absence of the interac- 
tion, J = 0, are very close, in average, to \ip n \ = l/\/N. 
If the interaction is very weak, the standard perturbation 
theory is valid, and a kind of fluctuations of the proba- 
bilities w n = \ip n \ 2 is expected around the mean value 
w n = l/N where N is the total size of the basis (the 
total number of many-particle states) . 

The data show that if the interaction J is relatively 
strong, the components of eigenstates are quite different 
from the unperturbed values. This region may be very 
important for quantum computation, and the main prob- 
lem is to know whether these errors in the components of 
the eigenf unctions (the deviations dw n from the unper- 
turbed value l/N) can destroy quantum coherent effects 
needed for the quantum computation. This problem was 
addressed in our previous study piij , here we are mainly 
interested in global properties of eigenstates for a very 
broad region of the interaction. 
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The most interesting conclusion which can be drawn 
from the numerical data for a weak enough interaction 
(see Fig.||(a-b)), is that the eigenstates turn out to have a 
regular structure, even if the deviations 5w n are relatively 
large. Indeed, one can see regular global dependence of 
w n on the basis number n, with some fluctuations around 
the mean. This fact seems to be directly related to the 
specific structure of the Hamiltonian matrix presented in 
Fig§ 




2048 4096 



2048 4096 



2048 4096 



ln> ln> ln> 

FIG. 8. Typical structure of eigenfunctions for differ- 
ent interaction strengths, J = 0, 0.0002, 0.1, 1, 10, 100 denoted 
by (a,b,c,d,e,f) respectively. Eigenstates are taken from the 
central energy band for L = 12, £1 = 100 , u>h — cjq + k, 
k = 1, L, luo = 100, v = wo, o = l. 



of an eigenstate in a given basis, is its "width" a(E) de- 
termined as 

1/2 



a(E) 



(23) 



Note, that in contrast to N pc which gives an effective 
number of large components, and is insensitive to the 
location of these components, the width a(E) does not 
"feel" the presence of "holes" in the sparse eigenstates. 
The latter fact can be used to distinguish chaotic ergodic 
states from the sparse ones. Namely, for fully extended 
but very sparse eigenstates, the value of a{E) is of the 
order of N, however, N pc is much less than N. 
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With an increase of the interaction, the regular struc- 
ture of eigenstates disappears and huge fluctuations of 
components of eigenstates emerge, see FigJ|(c-e). The 
structure of these eigenstates is very similar to that 
known in the physics of disordered systems, when eigen- 
states "occupy" some fraction of the basis, without no- 
ticeable correlations between different components ip n 
(see for example, [ p3| and references therein). One can 
say that these eigenstates are sparse in the sense that the 
number N pc of principal components of the eigenstates is 
much less than the total size N of the basis. Therefore, 
there is a strong change in the structure of eigenstates 
(compare Figs.^(a-b) with FigsJ|(c-e)). One can say that 
the transition from extended regular states to the weakly 
chaotic states occurs for J w 0.1. 

When the interaction between qubits increases further, 
one can see another transition to strongly chaotic states, 
see Fig.||(f). The latter are characterized by an ergodic 
filling of the whole basis, and by strong fluctuations of the 
components ip n , which are found to be practically ran- 
dom and independent. This situation is well described 
by RMT (see, e.g. fl32fl). Therefore, for such strong in- 
teraction J « 100, chaotic properties of our system are 
very strong and quantum computation process can be 
destroyed. 

In order to quantitatively characterize the eigenstates, 
we have computed the number N pc of principal compo- 
nents defined by Eq. ( |l7| ) . Another measure of the spread 



FIG. 9. Normalized average number of principal compo- 
nents N pc (open circles) and the width a (full circles) as a 
function of J in ^-representation for fi = 100. The average 
is taken over the eigenfunctions from the central band only. 
The solid horizontal line corresponds to N pc = N, and the 
dotted line gives the extreme limit of completely chaotic and 
extended states, N pc — N/3. The parameters are the same as 
in Figg 

The mean values of N pc and a in dependence on the 
interaction J are given in Fig.^. The circles represent 
the value of N pc and a, averaged over the eigenstates 
from the central energy band. First of all, one should 
note that the width a turns out to be large and indepen- 
dent of the interaction. This means that all eigenstates 
are extended in the z-representation, in spite of a seri- 
ous difference in their structure, see Fig.||. Contrary, the 
number of principal components N pc demonstrates two 
principal transitions in the structure of eigenstates. 

Numerical data of Figs.|||| allows one to distinguish be- 
tween few different regions of the interaction strength J. 
The first region with a very weak interaction J < 2 • 10~ 3 
is characterized by the constant value N pc w N and corre- 
sponds to completely extended (|VVi| 2 « ^/N) eigenstates 
shown in Fig.^(a-b). In this region the energy spectrum 
consists of many close quasi-degenerate levels, thus lead- 
ing to a strong deviation from the Poisson distribution, 
see Sect. III. 
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behavior of N pc close to the band edges, compare Fig.|10|d 
and Figjlp|e. Namely, in the region of parameters of 
Fig.|lO|d, the highest value of N pc corresponds to the band 
edges, in contrast to Figjl^e where at the band edges 
the eigenstates are extremely localized (with a very small 
value of N pc ). The origin of this difference is not clear, 
however, it should be noted that the data reported in 
Fig.|lO|c have already been observed (and explained) in 
few models of isolated systems with interacting particles 
(see, for example, Q and 34 1). For those models it 
was found that for the unperturbed eigenstates which 
are close to the band edges, the interaction with other 
basis states is strongly suppressed. 



FIG. 10. Number of principal components N pc for all 
eigenstates reordered in increasing energy (|0) is the ground 
state, |1) is the first excited state, etc.). Data correspond to 
the parameters of Fig.^. 

In the second region with N pc <C N, all eigenstates 
are strongly influenced by the inter-qubit interaction. 
This region was termed in Ref. |?8) the region of weak 
chaos since the structure of eigenstates looks chaotic (see 
Fig.^d), however, the level spacing distribution P(s) is 
quite close to the Poisson. From the data, the transition 
to the weak chaos occurs for J w 0.05 and corresponds to 
the analytical estimate ([l6]) for the derealization transi- 
tion in the mean-field basis. The very point is that the 
critical value J cr given byEq.((l|), in the z-representation 
corresponds to the transition from completely extended 
states to the weakly chaotic states. One should stress 
that from the practical point of view the region of weak 
chaos may be dangerous for quantum computation be- 
cause of large deviations of eigenstates from the unper- 
turbed ones, see Fig|8|(c-d). 

Second transition to strong quantum chaos occurs for 
J ~ 100. By the latter term we denote the situation when 
the level spacing distribution has the Wigner-Dyson form 
and fluctuations of components ip n are close to Gaussian 
ones with N pc w N/3, see Fig.|[ As we already discussed, 
this transition corresponds to the simultaneous occur- 
rence of both band overlapping and delocalized states, 
see (p2|). One can see that strong quantum chaos for N- 
interaction emerges for an extremely strong interaction 
and thus it is not relevant for quantum computation. 

More detailed information about the global structure 
of eigenstates can be drawn from Fig.|lQ where the value 
of Npc is shown for all eigenstates ip n (£^ m ) ) reordered in 
increasing energy E^ m '. In this figure one can see how 
the band structure of the spectrum manifests itself in 
the value of N pc . In particular, it is seen that for non- 
overlapped bands there is a quite strong dependence of 
N pc on whether the energy of a specific eigenstate is 
at the center of energy bands or close to the band edges. 

One should point out a remarkable difference for the 



VI. RANDOM MODELS 

In the previous Sections we have discussed the dynam- 
ical model (|J) of interacting qubits. We have seen that 
in spite of the absence of any randomness in this model, 
for a very strong interaction both energy spectra and 
structure of eigenstates reveal chaotic properties which 
are generic for quantum chaos. In this sense, it is in- 
teresting to compare the obtained results with those for 
similar models with random interaction. This problem is 
not academic since in reality there are many effects which 
can lead to some randomness in the Hamiltonian (||). 



A. All-to-all interaction 

It is instructive to see what happens for a long-range 
interaction between qubits. We have studied in details 
the case when the interaction couples all qubits in the 
same manner (A-interaction), 



L-l 



fc=0 



-s k F k + mi - 2 



TZ TZ 

1 k 1 n 



(24) 



Here the interaction is assumed to be completely random, 
with Jk.n = JS, where £ are random numbers with a flat 
distribution inside the interval [— 1, +1], 

This model can be treated analytically in the same way 
as we did it in Sect. IV. Specifically, we are interested 
in the derealization border which is determined by the 
comparison of the ratio ( |l4| ) with the typical interaction 
strength. 

The modification of the Hamiltonian (^) written in 
the mean field basis is straightforward. Specifically, the 
structure of the unperturbed part, see Eq.(ffl), remains 
the same, and the interaction term (^J) has the same 
structure (the only difference being the summation taken 
over all qubits). The most important point is that the 
Hamiltonian matrix has a different structure from that 
for the N— interaction, see FigJTl] 
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FIG. 11. Structure of the Hamiltonian matrix in the mean 
field basis for the A-interaction; black points stand for matrix 
elements whose modulus is larger that 10 -6 . Here is L — 8, 
n = 100, J = 1, U) k = U) + k. 

Despite the block structure shared by the analogous 
matrix for the A^-interaction, shown in Fig.(^|), and due 
to two-body interaction, each block is now characterized 
by many elements different from zero. For this reason, 
one can expect that chaotic properties of the model with 
the A-interaction are much stronger that those found in 
the case of TV-interaction. 

The estimate for Mf can be obtained for the A- 
interaction as well. Since all qubits are allowed to interact 
to each other, the maximum number of coupling between 
unperturbed many-body states inside the central energy 
band with all the others, is 

Mf = ^. (25) 




\o~' io° id id id id id 



J/J cr 

FIG. 12. The average number of principal components in 
the mean field basis, for the eigenstates from the central en- 
ergy band, as a function of J/J cr for L — 12 and f2 = 1000. 
Open circles are for the A-interaction, full circles are for the 
TViV-interaction, see next Section. 



This is an unexpected result since it coincides with the 
estimate ( |l6| ) for the derealization border in the case 
of A-interaction. The reason is that the energy range 
(AE) f within which many-body states are connected by 
the interaction, and the number Mf of the states within 
this energy range, are both proportional to L? . The re- 
sult shows that the derealization border turns out to be 
independent on the range of the interqubit interaction. 



As for (AE)f, it can be found by considering the 
maximal energy shift obtained by applying the opera- 
tor IqI v l _i to the state |1l-i, 0o), and resulting in the 
new state |0l-i, lo)- The energy difference between 
these two states is given by 



(AE )f * - { 2L% 

which perfectly agrees with the direct computations, see 
Fig|| As a result, the critical value J° r for the dereal- 
ization border is obtained from the relation, 



ja 
_ cr 



(Aff)/ _ 2a 2 



M 



therefore, 



J a 



4a? 

n 



(26) 



However, chaotic properties of this random model with 
the ^-interaction are much stronger than those found 
for the AMnteraction. Namely, the chaos border for the 
AMnteraction turns out to coincide with the derealiza- 
tion border. The transition to delocalized states for the 
A— interaction is shown in Fig.[l2]. 



The closeness of the derealization and chaos borders 
for the A-interaction can be also checked by studying the 
level spacing distribution. The latter is expected to man- 
ifest a transition from the Poisson to the Wigner-Dyson 
at the critical value of J given by the above estimate (prj). 
In Figjl^ we show that the transition to chaos is inde- 
pendent from the product Jf2, in correspondence with 
the analytical prediction (|2^). These results prove that 
for the A-interaction our model is similar to generic mod- 
els for which the derealization border coincides with the 
chaos border. 
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FIG. 13. Level spacing distribution for eigenvalues in the 
central band for L = 10 and 5k = k. For the average, 30 
different matrices with the random A-interaction have been 
used, a) f2 = 10, J = 0.1 (open circles), J = 1 (full cir- 
cles); b) D. = 100, J = 0.01 (open squares), J = 0.1 (full 
squares); c) Q = 1000, J = 0.001 (open diamond), J = 0.01 
(full diamond). Note that the theory predicts a transition 
point at J = J cr ~ 4/Q. For comparison, both the Poisson 
(dashed line) and the Wigner-Dyson (full line) distributions 
are shown. 



FIG. 14. Nearest neighbor distribution for eigenvalues in 
the central band for L — 12. One single matrix with ran- 
dom AW-interaction have been used. Open circles are for 
J — 0.001, close circles for J — 1. For comparison, both the 
Poisson (dashed line) and Wigner-Dyson (full line) distribu- 
tions have been shown. 

Since any other of long range interactions can be seen 
within these two extreme cases (the A and NN inter- 
actions), one can conclude that for a typical interaction 
(other than strictly between nearest qubits), the quan- 
tum chaos can emerge for a quite weak interaction and 
may have an influence on a quantum computer operabil- 
ity. Therefore, it may be important to reduce the range 
of the interqubit interaction in an experimental setup of 
a quantum computer. 



B. Next to nearest interaction 



VII. GENERAL DISCUSSION 



A. Quasi-integrability 



Finally, we discuss the intermediate case when the in- 
teraction V in the dynamical model (Q) couples four next 
nearest qubits, k ± 1, k ± 2, (the NN- interaction). 

A straightforward analysis similar to that shown in the 
previous Sections leads to the same critical border for 
delocalized states, as those found for the N and A inter- 
actions. This has been numerically confirmed, see that 
data in Fig.|l2|. Moreover, as for the ^-interaction, the 
delocalization border for the AA-interaction turns out 
to coincide with the chaos border. This has been proved 



by using the level spacing distribution, see Fig. 14 



Our numerical study shows that, in contrast to the 
case of the TV-interaction (when only two neighbor qubits 
are coupled), the quantum chaos emerges for much lower 
values of the AA-interaction, for 0.1 < J < 1.0, see 
FigsQ This region of parameters J and O is important 
from the experimental viewpoint, therefore, the quantum 
chaos may have a real influence for quantum computa- 
tion. 



As we noted, the model (|J) with the interaction 
between nearest qubits has quite specific properties. 
Namely, the delocalization border turns out to be very 
different from the border of quantum chaos. Below we 
explain this phenomena in terms of quasi-integrability of 
our model. 

Let us come back to the expression for the off-diagonal 
matrix elements of the Hamiltonian (||) in the mean-field 
basis determined by the eigenstates of H , see Eqs.(||). 
For the case of large O ^> 5k, we are interested in, the 
term V // is small compared to Vband since ak ~ 1 and 
bk ~ — I/O. Also, the diagonal term Vdiag is much 
smaller than the two other terms (it is proportional to 
b\ ~ I/O 2 <C 1/0 -C 1, see Eq.@). Therefore, the 
approximate Hamiltonian H a can be written in the fol- 
lowing form 



H a = 



L-l 

k=a 



L-2 



(27) 



k=0 



where 7fe = y/Sf. + O 2 and Jk = 2 J for our model. 
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This Hamiltonian has been recently studied in a num- 
ber of papers, (see, for example, |35|] and references 
therein) . It was shown [ ^6| that for independent random 
variables jk and the model (|27]) can be mapped to an 
Hamiltonian describing L free fcrmions. This transfor- 
mation holds only in the case of nearest neighbor cou- 
pling. Therefore, this model is integrable and the level 
spacing distribution P(s) can be expected to be Poisson- 
like for any interaction strength < J? > 1//2 . This explains 
why for non-overlapping bands our original Hamiltonian 
(^) with Q 3> Sk reveals the Poisson for P(s) above the 
derealization border. 

It should be noted that the derealization border J cr 
(see Sect. IV) results from the standard perturbation the- 
ory which takes into account a two-body nature of inter- 
action. Namely, when the typical interaction which con- 
nects unperturbed many-body states is much larger than 
the mean distance between energy levels of these states, 
in the corresponding basis the interaction creates ex- 
act eigenstates with many components. Typically, these 
compound states are chaotic due to a complex structure 
of the interaction. This is why the derealization bor- 
der generically coincides with the quantum chaos border. 
However, in specific cases like our quasi-integrable model 
(for Q 3> Sk and not very strong interaction), the de- 
localization border and the onset of chaos may be very 
different. 

The above analysis is also helpful in the explana- 
tion of the strong difference between the model with N- 
interaction, from the model when qubits are coupled by 
a different kind of interaction ( A— or NN— interaction, 
see previous Sections). Indeed, in the latter cases the 
interaction V has many additional terms compared with 
Eq.(||), and results in a strong coupling between all en- 
ergy bands. This leads to quasi-integrability breaking, 
and to the onset of chaos at the border of derealization. 



B. Role of magnetic field 



Our approach based on the mean-field representation, 
see Sect. IV, is valid for any kind of the B 2 -magnetic field. 
Let us consider the simplest case of a homogeneous mag- 
netic field for which all frequencies of spin's precession u>k 
are the same, StOk = u>o ~ v = /• For a non-resonant case 
with / ^ 0, and in absence of the interaction (J = 0), 
the energy spectrum has no more a band structure since 
each of the L+ 1 levels is degenerate. Indeed, each single- 
particle energy has two values e& = ±|(^ + ^n) only, 
where / < t]. Since all many-body states in the central 
band have the same number of pluses and minuses in the 
expression for the total energy, the latter is zero. Thus, 
the level spacing (AE) / is also zero which means that 
any small interaction gives rise to delocalized states. 



CO 2 

Z 



o A=0.1,L=8 
+ A=1,L=8 
x A=10, L=8 
♦ A=10,L=12 



♦♦♦ 



j®« ®«® See ©®ffi 6oo 

♦♦♦ -F 

x 

xxx^o<x 



Iff 
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FIG. 15. Average number of principal components for 
eigenfunctions in the central band for homogeneous magnetic 
field and random frequencies in the interval (f— A/2, v+A/2), 
versus the rescaled interaction J/J cr , where J cr is defined by 
Eq-dl). 

In recent studies p4j random variation of spin frequen- 
cies is included in the model, in order to take into account 
effects of finite temperature and environment. For this 
reason the energies are not exactly degenerate but swap 
into finite width bands. In the same way, let us assume 
that the energy of many-body states fluctuates, thus re- 
sulting in the distribution of the parameter / within some 
interval (-A/2, +A/2) with A < SI. Then, one can es- 
timate, 

A 2 

On the other side, the number of coupled state for a fixed 
state from the central band remains the same, Mf f=a L/2. 
As a result the derealization border can be determined 
from the relation, 



J > J r . 



(A£)y 
M f 



4QL 



(28) 



This parametric dependence has been checked numer- 
ically (see Fig.|l5|) where the average number of princi- 
pal component is plotted against the rescaled interaction 
J/ J cr for different L and A. As one can see, the scal- 
ing law given by Eq.(|28|) is quite well satisfied. Com- 
paring with Fig.||, one should note that for a constant 
magnetic field the onset of a strong chaos (N pc s» A/3) 
happens in a very small region of interaction (see the 
presence of small peaks on the far right side). With fur- 
ther increase of the interaction, the system again becomes 
nearly-integrable, since in the limit J 3> SI only diagonal 
terms dominate. 

In this way we come to the same L-dependence for 
the critical interaction J cr , discussed in Refs. In 
these papers, the model with a nearest interaction in the 
plane was considered (rather than on ID line as in our 
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model). For this reason the model of Refs. H| is free 
from the effects of quasi-integrability and, therefore, the 
derealization border coincides with the border of quan- 
tum chaos. 

Finally, we would like to point out that in the case of 
increasing gradient of the Z3 z -magnetic field, the dereal- 
ization border increases with an increase of the number 
of qubits. This very unexpected prediction can be easily 
understood for the case ojk = bk 2 (linear increase of the 
gradient). It can be shown that the width (AE)f grows 
proportional to L 3 , therefore, for the nearest interaction 
(Mf ~ L) the critical interaction increases as J cr ~ L 2 , 
and for the yl-interaction one gets, J cr ~ L. In the lat- 
ter case the estimate of J cr also gives the transition to 
the chaos. As one can see, the magnetic field with an 
increasing gradient may strongly reduce the influence of 
the derealization and chaos. 
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